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We consider inhomogeneous big bang nucleosynthesis in 
light of the present observational situation. Different obser- 
vations of 4 He and D disagree with each other, and depending 
on which set of observations one uses, the estimated primor- 
dial 4 He corresponds to a lower baryon density in standard 
big bang nucleosynthesis than what one gets from deuterium. 
Recent Kamiokande results rule out a favorite particle physics 
solution to this tension between 4 He and D. Inhomogeneous 
nucleosynthesis can alleviate this tension, but the more likely 
solution is systematics in the observations. The upper limit 
to Sit from inhomogeneous nucleosynthesis is higher than in 
standard nucleosynthesis, given that the distance scale of the 
inhomogeneity is near the optimal value, which maximizes ef- 
fects of neutron diffusion. Possible sources of baryon inhomo- 
geneity include the QCD and electroweak phase transitions. 
The distance scale of the inhomogeneities arising from the 
electroweak transition is too small for them to have a large ef- 
fect on nucleosynthesis, but the effect may still be larger than 
some of the other small corrections recently incorporated to 
SBBN codes. 



have reviewed the 4 He observations and their best esti- 
mate is 



PACS numbers: 



i.Cq, 98.80.Ft 



I. INTRODUCTION 



Standard big bang nucleosynthesis JT|-|| (SBBN) pre- 
dicts the primordial abundances of D, ^He, 4 He, and 7 Li 
as a function of a single parameter, the baryon-to-photon 
ratio rj = nf,/n 7 , which is related to the baryonic mass- 
density parameter Jlf, = 8irGpbo/3H 2 by 

n b h 2 = 3.70 x l(T 3 ?7io, (1) 

where 7710 = 10 10 r/ and h = ff /100kms _1 Mpc _1 . The 
observed abundances of these isotopes are in a rough 
agreement with the SBBN predictions || for a range of 
?7io, which is compatible with other cosmological bounds 
on the amount of baryonic matter in the universe. In 
principle, comparing SBBN predictions with primordial 
abundances extrapolated from observations pins down 
the precise value of ryig. A few years ago the standard 
result was 7710 ~ 3-4 |^|], but even much tighter con- 
straints were published (e.g., 2.69 < 7710 < 3.12 Q). Re- 
cently the situation has become more complicated, and it 
seems that such precise determinations were premature. 

Since the discovery of the r lepton, implying three 
flavours of light neutrinos, there has been tension be- 
tween 4 He and D in SBBN @J|. Olive et al. § (OSS97) 



Y p = 0.230 ±0.003. 



(2) 



This corresponds to 7710 = 1.4 ±0.3 and hence to primor- 
dial D/H ~ 2-3 x 10" 4 in SBBN, whereas the present 
D/H in the ISM is go) only 1.5 x 10" 5 . Most models 
of galactic chemical evolution have difficulty explaining 
this much deuterium astration [pd|, and prefer a much 
lower primordial D/H and thus a higher baryon density, 
7710 ~ 5. 

The conventional way to deal with this tension has 
been to compromise by settling on an intermediate 7710 
which is preferred neither by 4 He nor by D/H but is con- 
sidered acceptable to both. This however, leads to an ar- 
tificially high precision in the 7710 determination, because 
while the individual ranges in 7710 accepted by 4 He and 
D/H are wide, their overlap is narrow. Tension increased 
when data was subjected to more thorough formal sta- 
tistical analysis, culminating in a claim of a "crisis" in 
SBBN, by Steigman et al. JT2J , who concluded that given 
the existing data the overlap is in fact nonexistent. 

In the context of SBBN the resolution of this crisis 
requires either a revision of the picture of the galactic 
chemical evolution [l3| , so that much more deuterium 
astration can be accommodated jl4|], or a large system- 
atic error in the Y p determination from the observations 
P, p^|p6| . Indeed, based on a number of new 4 He observa- 
tions, Izotov and Thuan fTlj] have claimed a significantly 
higher Y p than OSS97: 



Y p = 0.244 ± 0.002. 



(3) 



Whether this new value is to be accepted as such over 
the old OSS97 value is yet unclear, since several sources 
of poorly known systematic effects are expected to con- 
tribute to the discrepancy p8| . 

Interestingly, some particle physics solutions based on 
a massive decaying tau neutrino Jl9| can now be ruled 
out using the recent results from Kamiokande pfj. The 
directional dependence in the upward going muon neu- 
trino deficiency seen in the Super Kamiokande experi- 
ment is a strong implication that the muon neutrinos 
undergo oscillations while traversing through the earth. 
This implies that mixes with either a tau neutrino 
or a new sterile neutrino with a mass splitting of about 
Sm 2 ~ 10~ 3 eV 2 and with an almost maximal mixing 
angle. If this mixing is between and v T , then v T is 
obviously light so that the scenarios based on heavy v r 
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decaying into and some scalar particle |Tst | are im- 
mediately ruled out. Suppose then that the atmospheric 
anomaly is due to mixing between za, and some sterile 
neutrino. Now v r can be heavy and having it decay away 
to muon neutrino and a scalar state prior nucleosynthe- 
sis could alleviate the tension somewhat. The effect is 
roughly equivalent to having about a half a neutrino de- 
gree of freedom worth less energy density in the universe 
Jl9[ (less energy density leads to slower expansion, and 
hence later decoupling of n/p-ratio). However, the ster- 
ile state with the requested mixing parameters is brought 
into full thermal equlibrium due to oscillation and quan- 
tum damping prior to nucleosynthesis pl[ , overcoming 
the alleviating effect discussed above and making the ten- 
sion even worse. The only possibility to alleviate the ten- 
sion is that m UT <, 1 MeV and v T decays into an electron 
neutrino in the short interval after the electron neutrino 
freezeout but prior the onset of nucleosynthesis. In this 
case the excess (almost thermal) electron neutrinos can 
significantly increase the weak interaction rates keeping 
the n/p-ratio in equilibrium longer and hence leading to 
much less helium being produced [ p2"|]23| ]. Bringing the 
sterile neutrino into equilibrium makes also this solution 
less effective, but is not strong enough to rule out the 
possibility entirely §|Q . 

The chemical evolution of D and 4 He is particularily 
simple: 4 He increases with time, whereas D decreases. 
In contradistinction, 3 He and 7 Li are both produced and 
destroyed during galactic chemical evolution. Thus it is 
much more difficult to make reliable claims of their pri- 
mordial abundances based on present abundances. Ob- 
served 3 He abundances |25| in particular vary a lot within 
the galaxy and 3 He observations are useful for constrain- 
ing BBN only when combined with D and chemical evo- 
lution models. 

For 7 Li there is a very impressive plateau [p6| of abun- 
dances in PopII stars with surface temperatures T — 
5800-6400 K. The observed value is 

log 10 ( 7 Li/H) = -9.75 ±0.10. (4) 

The universality of this abundance suggests that it is 
closely related to the primordial abundance. There may 
have been some depletion, i.e., some of the surface 7 Li 
has been destroyed by the star. Therefore the primordial 
abundance should be larger by some depletion factor D7. 
Pinsonneault et al. [§8[ estimate D7 = 0.2-0.4 dex. This 
corresponds to a primordial 

log 10 ( 7 Li/H) p = -9.45 ± 0.20. (5) 

However, Vauclair and Charbonnel |2{| give a lower esti- 
mate 

log 10 ( 7 Li/H) p = -9.65 ± 0.10. (6) 

These estimates for lithium are compatible with either a 
low, ijio ~ 1.5, or a high, ?/ 10 ~ 4-6 baryon density, but 
disfavor a compromise value rjio ~ 2.5-3. 



A promising new method with the potential to resolve 
this rj dicothomy is the observation of deuterium in clouds 
at high redshifts by its absorption of quasar light [ p0[ . 
Some of these clouds are so far away, that when the 
observed light passed through them, the universe was a 
mere one-tenth of its present age; thus the matter in these 
clouds and therefore the observed deuterium abundance 
should be close to primordial. Unfortunately, at present 
we only have a small number of such D/H measurements, 
and even the existing ones are still controversial. Buries 
and Tytler [ pl| obtain from their two best observations 

D/H = 3.4 ± 0.3 x 10~ 5 , (7) 

which corresponds to r/ia = 5.1 ± 0.3 in SBBN. This is in 
contradiction with the OSS97 estimate Y p ~ 0.23. How- 
ever, the analysis of Buries and Tytler has been debated 
p2| and one observation by HST |Q from an absorp- 
tion cloud at z = 0.7 appears to give a high value of 
D/H~2xl0~ 4 . 

Thus the observational situation remains unclear. If 
we suppose that some of the determinations of primordial 
abundances are correct, but we do not know which, we 
are led to an SBBN range 

1710 ~ 1-5-6. (8) 

One can also try to determine the universal baryon 
density by other means, discarding nucleosynthesis con- 
siderations |3^ ]. Determinations of this kind are rather 
uncertain at present, but tend to favor the larger values 

of 1]. 

In conclusion, there is an unsettled disagreement be- 
tween different observations in the context of SBBN. 
While the problem may lie with the observations, or in 
the determination of primordial abundances from them, 
another possibility is, that the primordial abundances in- 
deed do not correspond to the same 77 in SBBN, so that 
it needs to be modified. In this paper we study the possi- 
bility of inhomogenous big bang nucleosynthesis (IBBN) 
in light of the present observational situation. In section 
II we discuss the generic mechanisms known to produce 
inhomogeneities in the baryon distribution and the sig- 
nificance of the distance scale of the inhomogeneity. We 
describe our numerical calculations in section III and give 
our results in section IV. Section V contains our conclu- 
sions. 



II. INHOMOGENEOUS NUCLEOSYNTHESIS 

In SBBN we assume baryonic matter was homoge- 
neously distributed during nucleosynthesis, but actually 
we do not know whether this was the case. If the inho- 
mogeneity was of sufficiently small scale, diffusion would 
have homogenized the matter distribution before the for- 
mation of the cosmic microwave background leaving no 
directly observable trace today. 
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A. Generating the inhomogeneity 

Various phase transitions which took place before nu- 
cleosynthesis in the early universe were capable of pro- 
ducing large-amplitude small-scale fluctuations in the 
baryon number density: in particular the electroweak 
(EW) transition at T ~ 100 GeV and t ~ 10~ n s and 
the QCD transition at T ~ 150 MeV and t ~ 10" 5 s. 

IBBN was studied extensively in the late 1980's, when 
it was realized that a first-order QCD phase transition in 
the early universe could produce the kind of inhomogene- 
ity which would affect BBN |^@. At first ||,||,|^ 
it seemed possible to accommodate much larger values 
of 7710 , even fij, = 1, but more detailed calculations 
|38|,fl0 42 4^ [ItJ showed that the upper bound to 7710 was 



in fact much less increased. 

The original mechanism relying on chemical pressure 
psfl , operative in the QCD transition, leads to a geometry 
where localized clumps of high density are surrounded 
by large voids of low baryon density p9|-p2| . The details 
of the QCD transition are poorly known and both the 
amplitude and the size of the inhomogeneities can vary 
significantly; the size of course is bounded by the horizon 
at the QCD transition, which is about 2 x 10 6 m (at T = 
1 MeV) = 0.4 pc (today). 

Also the electroweak phase transition (EWPT) gener- 
ically produces inhomogeneities, and possibly with large 
density contrasts. This assumes of course that the 
baryons we see around us today, were generated dur- 
ing the electroweak phase transition [Q. Some scenar- 
ios [p4 55 1 may even give rise to regions of antibaryons 
mixed with the overall baryonic excess, leading to the in- 
teresting possibility of nucleosynthesis in the presence of 
antibaryons [ p6|j57| . The generic feature leading to the 
formation of inhomogeneities in the more standard sce- 
narios is the strong dependence of the baryoproduction 
rate on the bubble wall velocity in the so called "charge 
transport mechanism" (5^-Q , coupled with the charac- 
teristic changes in the velocity of the bubble walls during 
the transition |6l[]. For thin walls one finds a local bary- 
oproduction rate 



B(x) ps c/v w (x). 



(9) 



The velocity dependence of the local baryoproduction 
rate due to the " classical chiral force" mechanism ^||,^3 1 , 
operative in the limit of wide walls, is much weaker [p4 1 . 
However, the generic geometry of inhomogeneities aris- 
ing from EWPT is quite opposite to the QCD case; voids 
of low density surrounded by walls of high density. 

After nucleation bubble walls quickly accelerate to a 
terminal velocity v w ~ 0.1-0. 5c, whose exact value de- 
pends on the parameters of the phase transition, like 
the latent heat released, the surface tension and the fric- 
tional forces effected on the bubble wall by the ambient 
plasma |65|,|66|]. After some time (we are only consider- 
ing deflagration bubbles here), the shock waves preced- 
ing phase transition fronts collide reheating the unbroken 



phase plasma back to the critical temperature. As a re- 
sult the pressure forces driving the bubble expansion are 
reduced, and, were it not for the general expansion of 
the universe, the walls would become to a complete stop. 
Due to Hubble expansion the walls can still continue ex- 
panding, but now with a greatly reduced speed, typically 
v w ~ O(few) x 10~ 3 c |66|. These velocity scales and the 
rate (||) indicate that the maximal density contrast pos- 
sibly generated by the EW mechanism is about ~ 100. 

The typical size of the voids in this "beer foam" ge- 
ometry is some fraction of the horizon at the EW transi- 
tion, £ H = 3 x 10 3 m (at 1 MeV) = 6 x 10~ 4 pc (today). 
A nucleation calculation, which ignores the thermody- 
namics of the bubble interactions, typically gives for the 
size of bubbles at the coalescence only If, ~ 10~ 3 £h 
p(i| |G8[ . However, due to reheating the firstly nucleated 
bubbles may inhibit the growth of bubbles formed only 
slightly later, increasing perhaps significantly the size of 
the largest structures as compared with the simplest nu- 
cleation estimate. Also in extended scenarios including 
magnetic fields 55 1, the size of a single bubble can reach 
the horizon scale. We then consider the inhomogeneity 
size a free parameter, with values r ~ 10~ 3 — 1£h- 



B. Distance scales 

Both the EW and the QCD transition appear capable 
of producing high initial density contrasts. In both cases 
the density fluctuations would be non-gaussian, consist- 
ing of high- and low-density regions. The pattern would 
not be regular, but it would have a characteristic distance 
scale. The inhomogeneity can be described by the typical 
geometric shape of these regions and the following three 
parameters: 1) typical distance scale r, 2) typical density 
contrast R = ?7high/^iow, and 3) the volume fraction /„ 
of the high-density regions. 

The distance scale r is especially important. An 
inhomogeneity can have a large effect on nucleosyn- 
thesis only if the distance scale is comparable to the 
neutron diffusion length d n during nucleosynthesis. If 
the distance scale is too small, r <§; <i„(500 keV) ~ 
200 m (at 1 MeV) ~4x 10" 5 pc (today), the inhomo- 
geneity is erased before nucleosynthesis, because before 
the weak freeze-out protons and neutrons are constantly 
converted to each other by weak reactions, and the dif- 
fusion thus evens out both the proton and the neutron 
density. 

If the distance scale is large, r ^> d n (10 keV) ~ 
500 km (at 1 MeV) ~ 0.1 pc (today), diffusion does not 
occur until nucleosynthesis is completed. In this "ordi- 
nary inhomogeneity" scenario, the high- and low-density 
regions undergo independent standard BBN with /Thigh 
and ?7i ow and the matter is mixed afterwards to have the 
average baryon-to-photon ratio 77. Leonard and Scherrer 
|69| have shown that this kind of inhomogeneity cannot 
increase the upper limit to 77, since the inhomogeneity 
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raises 4 He yields. Arbitrarily low 77 can be made accept- 
able with these models, however. 

An intermediate distance scale, d„(500 keV) < r < 
d n (10 keV) leads to nucleosynthesis with inhomogeneous 
neutron-to-proton ratio. The strongest effect occurs if 
n/p 3> 1 in some region, because in this region neutrons 
are left over from 4 He synthesis. This can be induced 
in the low-density region if f v R ^> (n/p)^ 1 before dif- 
fusion. QCD-scale inhomogeneities could be [plj of the 
scale required, although QCD lattice calculations |7(| fa- 
vor values below the short end of this range. For the 
EW case this range corresponds to a fluctuation scale 
r ^ 0.1£h during the transition. 

There may be other possible sources of baryon inho- 
mogeneity in addition to the EW and QCD phase tran- 
sitions. Moreover, there is a considerable uncertainty 
regarding the parameters r, R and from each transi- 
tion. Therefore it is natural to treat the two questions 
separately: 1) Are there IBBN parameter regions where 
IBBN agrees with observations equally well or better 
than SBBN? 2) Could the EW or QCD phase transitions 
produce inhomogeneity in this parameter region? 



III. COMPUTATIONS 

The IBBN code used for this paper is based on the 
code used in jl3| and the nuclear reaction rates have been 
updated according to In the 4 He yield we take into 
account the various corrections to the weak reaction rates 
P,p|j7l|] . Theoretical uncertainty in abundance yields due 
to uncertainty in nuclear reaction rates is usually small 
compared to observational uncertainties. An exception 
is the 7 Li yield for which one standard deviation is 0.07- 
0.10 dex upwards and 0.11-0.19 downwards 0. We take 
this into account when obtaining limits on 77 from 7 Li 
deduced from observations, by further relaxing the up- 
per limit to 7 Li by 0.15 dex and the lower limit by 0.10 
dex. Proton diffusion is included according to |7^| . The 
convergence of the code has been improved by combin- 
ing the nuclear reaction and diffusion steps into a single 
step, which reduces the number of time steps needed for 
accurate results. 

We assume spherical symmetry and use a nonuniform 
radial grid of 64 zones representing a sphere with comov- 
ing radius r, with reflective boundary conditions both 
at the center and at r. This setup allows us to model 
both geometries discussed above: assuming centrally con- 
densed density describes the QCD-type and spherical 
shells of high density describe the EW-type geometry. 
The volume fraction covered by the high density region 
in each geometry is 



fv f r 

fv = l-(l- fr? 



(centrally condensed), (10) 
(spherical shell), (11) 



is specified by four parameters: r, /„, R described above, 
and the average baryon-to-photon ratio 77. The numerical 
value for 77 always refers to the present value, i.e., after 
e + e~-annihilation. 

Note that these inhomogeneities are in baryon num- 
ber only. At nucleosynthesis time, the energy density is 
dominated by radiation by a factor of at least 10 5 . The 
density of baryon number can hence be strongly inho- 
mogeneous without a noticeable dynamic effect, and the 
main process through which the inhomogeneity evolves 
after the phase transition is diffusion [^,Q|6). 

The effect of neutron diffusion is to reduce 4 He and 
7 Li yields, and to increase D and 3 He yields. All these 
changes are in the direction of favoring a larger 77. How- 
ever, diffusion has to compete with the ordinary inhomo- 
geneity effect which for 4 He, 7 Li, and 3 He is the opposite, 
increasing 4 He and 7 Li and reducing 3 He. For D this lat- 
ter effect depends on the average 77. For small 77, D is 
reduced, and for large 77 D is increased. 

The most dramatic effect is obtained when the neu- 
tron diffusion out of the high-density region leads to a 
large excess of neutrons in the low-density region. This 
requires a density contrast 



P) 1 



(12) 



where (p/n)o ~ 7 is the SBBN proton/neutron ratio at 
the onset of nucleosynthesis. Increasing R leads to a 
stronger effect, but the increase soon saturates. Indeed, 
for large R almost all of the nuclear matter already is in 
the high-density region, while almost no matter remains 
in the low density regions: 



77high = 
'/low 



Rr] R-^oa J]_ 

f v R + 1 — /„ f v 

V 77 ft- 

f v R+l-f v ~ JJi 



0. 



(13) 
(14) 



R 


fr 


fv 


fvR 


Centrally condensed (c.c.) 


283 


1/V2 


0.3536 


100 


800 


1/2 


0.125 


100 


2263 


1/2^2 


0.0442 


100 


6400 


1/4 


0.0156 


100 


51200 


1/8 


0.0020 


100 




Spherical shell (1 






1000 


1/4 


0.5781 


578 


1000 


1/8 


0.3301 


330 


1000 


1/16 


0.1760 


176 


1000 


1/32 


0.0909 


91 


1000 


1/64 


0.0461 


46 


1000 


1/128 


0.0233 


23 



where f r denotes the fraction of the radius covered by 
the high-density region. Given the geometry, the model 



TABLE I. The different geometries studied. R is the den- 
sity contrast between the high and low density, f r is the 
high-density fraction of the grid radius, and f v is the cor- 
responding volume fraction. 
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The effect of further increasing R, beyond, say R = 
100/ /„, just leads to a further reduction of matter density 
in the low-density region and has essentially no effect on 
nuclear yields. An exception to this may be D, since the 
D yield drops so fast with increasing 77, that a significant 
part of the D yield may still come from the low-density 
region, giving rise to sensitivity to a reduction of 7710W and 
hence to R. 

In most cases we chose to run with large enough R to 
have close to this maximal effect. This leaves us with 
three parameters /„, 77, r. We did runs with 11 different 
values of /„ altogether (Table I). 

For the runs with spherical shell geometry, we kept 
R = 1000 constant. For the centrally condensed geom- 
etry some of the volume fractions were so small, that a 
larger R was needed to get the large inhomogeneity ef- 
fect. For the centrally condensed runs we kept the prod- 
uct f v R = 100 constant instead. 

IV. RESULTS 

It has been customary in IBBN studies 
to plot the regions in the (77, r)-plane allowed by different 
observational constraints. Since the observational situa- 
tion has become rather less clear recently, we present the 
results first as abundance contours for a given geometry 
and /„ , so different constraints can then be applied after- 
wards. For 4 He we plot the mass fraction Y p , for D and 
7 Li we plot the number ratios D/H and 7 Li/H. To save 
space, the less interesting 3 He is not shown. 

In Fig. [l] we show the results for the spherical shell 
(s.s.) geometry with R = 1000 and f r = 1/16. It is 
clear that at the distance scales attainable in the EW 
transition (indicated by the lower horizontal dashed line 
in the figure) the IBBN results do not significantly differ 
from SBBN results; the observational uncertainties are 
certainly much larger. However, even with scales as small 
as r r-j 0.05£h, the effect of inhomogeneity (see Fig. ||) can 
be larger than certain small corrections recently included 
into the SBBN computations | [fi| . 

Now we take a somewhat different point of view and 
consider a broader range of density contrasts and distance 
scales than can be produced in the electroweak phase 
transition. In Figs. ]| and ^ we present the results for 
the centrally condensed (c.c.) runs with f r = 1/2 and 
f r = 1/4. The results from other runs described in Table 
I are qualitatively similar. 




FIG. 1. The 4 He, D, and 7 Li yields from inhomogeneous 
nucleosynthesis runs with the spherical shell geometry, with 
R = 1000 and f r = 1/16 (/„ = 0.176). The contours of (a) Y p , 
(b) log 10 D/H, and (c) log 10 7 Li/H are plotted as a function of 
the average baryon-to-photon ratio 77 and the distance scale r 
of the inhomogeneity. The two horizontal dashed lines denote 
the horizon scale £h at the QCD (upper) and EW (lower) 
phase transitions. 
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FIG. 2. Effects of small scale inhomogeneity on the isotope 
yields. This figure is for the spherical shell geometry appro- 
priate for the EW transition, and for f v = 0.3301 and r)io = 5. 
The three lines correspond to R — 10 (solid), 100 (dashed), 
and 1000 (dotted) (or f v R — 3.3,33,333), showing how the 
effect saturates for large R, so that there is little difference 
between R — 100 and R — 1000. The horizontal axis gives 
the ratio of the distance scale to the EWPT horizon. 



A. Optimum scales 

How much 4 He is produced depends on the number of 
neutrons available. The yield is minimized at an opti- 
mum distance scale r opt ~ 10 4 — 10 5 m, where a maximal 
number of neutrons diffuse out from the high-density re- 
gion (where most of the 4 He is produced), but not too 
many of them diffuse back when the nucleosynthesis in 
the high-density region starts consuming free neutrons, 
and the direction of the neutron diffusion reverses. D 



yields are maximized at scales somewhat larger than r op t , 
in particular for large 77, because of the strong ordinary 
inhomogeneity effect. Situation is more complicated with 
the 7 Li yields, but they tend to be minimized at r <, r opt . 




FIG. 3. Same as Fig. |l| but for the centrally condensed 
geometry with R = 800 and f r = 1/2 (/„ = 0.125). 
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FIG. 4. Same as Fig. but for the centrally condensed 
geometry with R = 6400 and f r = 1/4 (/„ = 0.0156). 

We find that r opt goes down with increasing 77, roughly 
as ?i~ 2 / 3 . Also, the s.s. geometry gives a larger optimum 
scale than does the c.c. geometry with the same f v , and 
the dependence on /„ is different with different geome- 
tries: for centrally condensed spheres r op t goes down with 



decreasing /„, whereas for spherical shells it increases 
with decreasing f v . 

It is possible to derive the parametric dependence of 
r opt on i] and /„ analytically. Consider the diffusion of 
neutrons after the weak freeze-out but before the start 
of nucleosynthesis. The flux of neutrons into the low- 
density region is proportional to the neutron diffusion 
coefficient D, to the surface area A of the boundary, and 
to the gradient of the neutron density at the boundary, 
roughly (rihigh — ni ow ) / V Dt. Here \f~Dt is the diffusion 
length of neutrons, and rihigh and n\ ow are the average 
neutron densities in the high- and low-density regions, 
respectively, rihigh decreases as 



V 



dn high 
dt 



AD 



T^high ^-low 



(15) 



where V is the volume of the high-density region. If we 
ignore nuclear reactions and weak inteactions, we can in- 
tegrate out Eq. (^5|). Remembering that /„ rihigh + (1 — 
fv)n\ ow = n mean one readily finds that the density con- 
trast vanishes exponentially due to diffusion: 



rihigh - "low ~ exp 



.4 



V (1 — f v ) 



The optimum scale corresponds to 



a yptz 

V(l- f v ) 



(16) 



(17) 



where t ns is the starting time of nucleosynthesis. At 
scales larger than the optimum scale, the neutrons have 
not diffused out effectively before the synthesis of 4 He 
begins. On the other hand, making the scale smaller 
than the optimum scale does not significantly increase the 
number of neutrons diffusing out, but makes the back- 
diffusion at later times more effective. 

Now it is easy to see why the optimum scales are 
smaller for condensed spheres. For the same /„, the 
surface-to- volume ratio A/V is smaller for condensed 
spheres than for shell geometry, which makes the out- 
diffusion less effective and optimum scales smaller. 

The efficiency of the out-diffusion also depends signif- 
icantly on the volume of the low-density region (term 
(1 — f v ) in the denominator of Eq. (p^)). If f v is 
large, n\ ow increases rapidly, bringing the diffusion to end 
sooner than in the case of small f v . 

The ry-dependence of the optimum scale is through the 
dependence on the diffusion length. The diffusion at the 
boundary is controlled by the smaller diffusion coefficient 
of the high-density region. The diffusion is dominated 
by scattering on protons, D np < D ne . After electron- 
positron annihilation the diffusion constant depends on 
the proton density and temperature as 



D r , 



D r< 



1 



oc 



fv 



hT 5/2 vT 5/2 



(18) 
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The early universe expands as t oc 1/T 2 . The diffusion 
length at the beginning of nucleosynthesis should go as 



V Dt ns OC %ji Tj 



-l/2 T _ 9 /4 



(19) 



The starting temperature T ns of nucleosynthesis depends 
on 77hi g h: in higher density nucleosynthesis begins earlier. 
The dependence in the range ?7high = 10~ 10 ---10~ 8 is 
T ns oc rg igh , 7 = 0.07 • ■ ■ 0.1. The diffusion length should 
depend on ??high as 



\J Dt us oc n^-li 2 97/4 ii, 



(20) 



l high — "high' 

where a = 0.65-0.73 ~ 2/3, so that for the optimum 
scale 



V 



(1-/^)2 «%igh 3 - 



(21) 



The surface-to- volume ratio of the high-density region is 
A 3 



A = 3 (1 - /„ ) 
V r f v 



(c.c) 



(22) 



2/3 



(s.s.). 



Combining equations ( pl| ) and ( |22| ) we find the observed 
behaviour for the optimum scale 



,1/3 

^o P t oc (c.c.) 



(23) 



^opt OC 



-2/3 



/, 1/3 (i"/,) 1/3 



B. Effects on the different isotopes 

4 He. For distance scales near the optimal one 4 He 
yields are reduced. The range in r where the He yield is 
below the SBBN value covers 1-1.5 orders of magnitude. 
For the optimum distance scale a given 4 He yield is ob- 
tained for a value of 77 which can be as much as four times 
larger than in SBBN. Making f v smaller causes a deeper 
reduction in He. For the spherical shell geometry we get 
the maximum effect at /,. ~ 1/64 although this may be 
due to our keeping R fixed to R = 1000. In the centrally 
condensed runs where we keep f v R constant instead, the 
effect keeps getting stronger for smaller /„. 

D. The runs near the optimum scale produce more D 
than SBBN. At these scales the D contours are pushed 
towards larger 77 by about the same amount as the 4 He 
contours. The maximum effect of the inhomogeneity is 
however at larger scales. Both D and 4 He consistently 
allow larger 77 in IBBN near the optimum distance scales. 

7 Li. Neutron diffusion helps to reduce the 7 Li yield, 
since at late times neutrons are diffusing into the high- 
density regions and destroying 7 Be there. However, the 
ordinary inhomogeneity effect is to increase 7 Li yields, 
and usually this effect wins, but for some cases we get a 
net reduction. 



C. Constraints on 77 

We now compare our IBBN yields to observational con- 
straints. Since at present there is no agreement about 
what constraints to use, we consider a number of differ- 
ent sets of constraints. 

The most fundamental abundance constraints are the 
upper limit to primordial He and the lower limit to 
primordial D/H, obtained directly from observed abun- 
dances, since chemical evolution always increases the 4 He 
abundance and reduces D/H. So, in our first set we con- 
servatively take for 4 He the 2-cr upper limit by Izotov and 
Thuan @, 

Y p < 0.248 (24) 

and for D/H we use the present ISM abundance jl0| as 
the lower limit, 

D/H > 1.5 x 10~ 5 . (25) 

It turns out that all our IBBN models which satisfy 
Eq. ( p3| ) satisfy Eq. ( |25| ) also. In Fig. || we have plotted 
the contour (£4]) from all the models considered here. 
In SBBN the constraint (|24|) gives an upper limit to 77, 
Vw < 6.3. We see that, e.g., the centrally condensed 
IBBN models with f r = 1/8 raise this upper limit to 



Vio < 19, 



or n b h 2 < 0.07. 



(26) 



Similar results were obtained for the spherical shell ge- 
ometry with f r = 1/32 or 1/64; reaching this upper limit 
requires the distance scale to be close to optimal in order 
to maximize the effect on rj. 

While IBBN raises the upper limit to 77 from 4 He and 
D/H by a factor of 2 to 3, upper limits from 7 Li are 
raised at most by a factor of 1.4, and, if we choose a 
very tight 7 Li limit, not at all. Thomas et al. used 
7 Li/H < 1.4 x 10~ 10 , which gives them an SBBN upper 
limit 7710 < 3.1, and this limit was not relaxed at all by 
IBBN. We confirm that none of our IBBN models raises 
the upper limit to 77 from this constraint. However, their 
upper limit for 7 Li allows essentially no depletion at all. 

As our second set we take the case for a high rj based 
on the high-z deuterium value of Buries and Tytler |jl| . 
We use the 2-cr range 



D/H = 3.4 ± 0.6 x 10" 



(27) 



as our constraint. For He we continue to use the Izotov, 
Thuan |l^] upper limit Y p < 0.248 and for 7 Li we use the 
Pinsonneault et al. |2q] range 



log 10 ( 7 Li/H) p = -9.45 ±0.20, 



(28) 



further relaxed by the theoretical uncertainties as dis- 
cussed in Sec. Ill, The results for this set are displayed 
in Fig. |. 
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In SBBN these constraints lead to a baryon density in 
the narrow range 7710 = 4.6-5.8. In IBBN the allowed 
range is 

??io = 3.9-8.2 (29) 
for the centrally condensed geometry and 

?/io = 3.7-10.5 (30) 
for the spherical shell geometry. 

u2 




10 t] 

FIG. 5. Conservative upper limit to r\ from Y p < 0.248 and 
D/H > 1.5 x 10 -5 . The plot (a) is for the c.c geometry: the 
thick curves are for f r = l/y/2 (solid), 1/2 (dashed), 1/2^2 
(dot-dashed), and the thin curves are for f r = 1/4 (solid) and 
f r = 1/8 (dashed). The plot (b) is for the s.s. geometry: 
the thick curves are for f r = 1/4 (solid), 1/8 (dashed), 1/16 
(dot-dashed), and the thin curves are for f r = 1/32 (solid), 
1/64 (dashed), and 1/128 (dot-dashed). The allowed region 
is to the left of each curve. 
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FIG. 6. The regions in the (r, r?)-plane allowed by D /H = 
3.4 ±0.6 x 10" 5 , Y p < 0.248, and log 10 7 Li/H = -9.45 ±0.20. 
Plot (a) is for the c.c. geometry and plot (b) is for the s.s. 
geometry. The meaning of the different line styles is the same 
as in Fig. H. 



In our third set we consider the case for low r\ in SBBN 
|§. (See Fig. 0). The 2-cr OSS97 limits for Y p 



0.224 < Y p < 0.236, 



(31) 



correspond to 7 Li near the Spite plateau and a large pri- 
mordial D. Hence we here use a conservative upper limit 
to D, 

D/H < 2.5 x 10~ 4 , (32) 
and the Vauclair, Charbonnel ]2Sf] upper limit for 7 Li, 
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log 10 7 Li/H < -9.55. 



(33) 



The results for this set are given in Fig. 0. The SBBN 
range is 7710 = 1.5-2.1 (lower limit from D/H, upper limit 
from Y p ). The IBBN upper limits are higher: 



Y p < 0.238, 
D/H < 10~ 4 , 



(SBBN 7710 < 2.4) 
(SBBN 7710 > 2.6) 



7710 < 3.6 (c.c, f r = 1/2) 
7/10 < 3.8 (s.s., f r = 1/16). 



(34) 
(35) 




10 ^ 10 

FIG. 7. The case for low 77. This figure is similar to 
Fig. [], but the constraints used are Y p — 0.230 ± 0.006, 
D/H < 2.5 x 10~ 4 , and log 10 7 Li/H < -9.55. 



We finally demonstrate that IBBN can alleviate the 
tension between low 4 He and low D. If we use the con- 
straints 



log 10 Li/H < -9.25, (SBBN r] W < 8.3) 



(36) 
(37) 
(38) 



no value of 77 is allowed in SBBN (the "crisis"). However, 
as is shown in Fig. ||, some IBBN models satisfy these 
constraints, with 2.6 < 7710 < 6.0 (c.c.) or 2.3 < 7710 < 6.5 
(s.s.), in a narrow (about a factor of two) range of the 
inhomogenity distance scale r. This is the "optimum" 
distance scale, which for these values of 77 varies between 
5 km and 30 km (at 1 MeV) for the centrally condensed 
geometry. Similar solutions were found with r about 70 
km for the spherical shell geometry, proving that the re- 
sult essentially depends only on the scale and is robust 
against using different geometries. 




10 ^ 10 

FIG. 8. Alleviating the BBN "crisis". This figure is similar 
to Figs. [] and but the constraints used are Y p < 0.238, 
D/H < 10~ 4 and 7 Li/H < 10~ 9 ' 25 . These constraints are 
incompatible with each other in SBBN, but are compatible in 
IBBN with the optimal distance scale. The allowed regions 
for the two geometries are shown in the same plot. The ones 
for the s.s. geometry are all at the same distance scale and 
lie on top of each other. For the c.c. geometry we get allowed 
regions for three of the considered volume fractions, and they 
lie below the s.s. regions. 



V. CONCLUSIONS 

In conclusion, we have studied the possibility of inho- 
mogenous nucleosynthesis on the basis of the new ob- 
servational situation, paying attention to the particular 
mechanisms capable of producing the inhomogeneities in 
the very early universe. 

First we studied the typical foam like inhomogeneity 
generated during the electroweak phase transition, which 
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we modelled by using spherical symmetry with thin shells 
of high density regions. We find that the scale from the 
EW transition tends to be too small to cause large devi- 
ations from SBBN predictions; that is, the bound on rj is 
not significantly changed. However, the effects on theo- 
retical yields can be of equal size or larger than some of 
the more detailed corrections recently included into the 
SBBN computations. Due to the genericity of the EW- 
inhomogeneities these corrections can be claimed to set 
the scale of accuracy achievable in SBBN computations. 

Second we considered the full parameter space of the 
IBBN models in both centrally condenced (QCD-type) 
and spherical shell (EW-type) geometries. 

To answer the first question posed at the end of Sec. |J: 
IBBN models can satisfy the observational constraints 
equally well, and for some small region of the parameter 
space, even better than SBBN. For inhomogeneities with 
distance scales near the "optimum" scale r op t , where the 
inhomogeneity effects are maximized, this agreement is 
obtained for a larger baryon density than in SBBN; pre- 
cise values depend intrinsically on the observational con- 
straints, but the upper limit to r/ from the upper limit to 
4 He and from the lower limit to D/H may be raised by a 
factor of 2-3, whereas upper limits set by the 7 Li/H data 
are raised less, at most by a factor of 1.4. However, it is 
not possible ever to get rj large enough to make Qt, = 1. 
For smaller scales the agreement is obtained for similar 
or slightly smaller values of 77 as in SBBN. 

Regarding the second question, this optimum distance 
scale is not only larger than the EWPT horizon, but it 
is also several orders of magnitude larger than the QCD 
transition distance scale favoured by QCD lattice calcula- 
tions of the surface tension and the latent heat. However, 
the uncertainty in these values is as large as the values 
themselves so that a much smaller latent heat, leading 
to a larger distance scale, is still allowed; thus we cannot 
presently rule out the possibility of reaching the optimum 
inhomogeneity distance scale in the QCD transition. 

There is a region of parameter space, where the tension 
between 4 He and D/H is alleviated compared to SBBN. 
This takes place if the inhomogeneity distance scale is 
close to r opt . The effect is however rather small, and 
for a low deuterium, say D/H < 5 x 10 -5 , we cannot ac- 
commodate less helium than Y p = 0.240, so IBBN cannot 
present itself as a solution to a dicothomy in observations. 
Since we also pointed out that the present Kamiokandc 
result rules out the simplest particle physics solution to 
possible tension in SBBN, the conclusion, that the prob- 
lems are probably associated with the observations, is 
bolstered. 
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